
(NASA-TM-X-65908) STABILITY OF 
MAGNETIC STARS AND PLANETS AS J 
THE ORIGIN OF CONVECTIVE MOTION 
(NASA) Jun. 1972 13 p 


V 1 ' 







STABILITY OF ROTATING MAGNETIC STARS AND PLANETS 
AS JUPITER AND THE ORIGIN 
OF CONVECTIVE MOTION 

Kunitomo Sakurai* 

Radio Astronomy Branch 
Laboratory for Extraterrestrial Physics 
NASA/Goddard Space Flight Center 
Greenbelt, Maryland, U.S.A. 

ABSTRACT 

Based on the method of the energy principle, the effect 
of the Coriolis force on the stability of rotating hydromagnetic 
systems is examined and the condition for instability is 
derived. It is shown that, in the rotating systems, the 
effect of this force is to inhibit the onset of convective 
motion and that, once the condition for instability is ful- 
filled, over-stable states are produced with respect to 
this motion. Such states seem to be effective on the onset 
of turbulent convective motion in the systems. 
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X. INTRODUCTION 

When a rotating magnetic star or planet is in hydrostatic 
equilibrium, this state is expressed by 

pfi x (fi x r) = - V P + — j x B + pg (1) 

C 

curl B = ~ j (2) 

and 

div B = 0, (3) 

where p, 0, r, P, c, j, B and g are the mass density, the 
angular velocity, the position vector, the pressure, the 
speed of light, the electric current, the magnetic field 
and the gravity force, respectively. 

If we take into account the finiteness of the electrical 
conductivity of the system, though very high, it is shown 
that the magnetic fields inside the system would necessarily 
decay gradually as a result of the Joule dissipation. The 
time for this decay is given by 

4 it o 2 . 

t = — £ L , (c.g.s.e.s.u) 

c 

where a and L are the electrical conductivity and the 

characteristic length of the system under consideration, 

respectively (e.g., Cowling, 1946, 1953). If we assume that 

12 

the electrical conductivity c is numerically given as ~10 
e.s.u., for instance, in the interior of the earth and Jupiter, 
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the above times for these two planets are estimated as 

4 0 i 

~10 and ~10 years, respectively. These magnitudes are 
too short to explain the magnetism of the earth and Jupiter 
by using the idea of "fossil" magnetism as proposed by 
Cowling (1946) . 

In order to maintain the magnetism of the earth, for 
example, we should, therefore, consider the re-generation 
process of magnetic fields in its interior. As we have 
shown above, it is evident that the system, in which there 
is no fluid motion, is not capable of preventing its own 
magnetic field from decaying completely. Some type of 
convective motions, which seem to be effective on re- 
generation or amplification of magnetic fields, must, there- 
fore, exist in order to explain the origin and maintenance 
of the magnetic fields of stellar objects such as the sun, 
Jupiter and the earth. 

2. S TABILITY IN THE NON-ROTATING SYSTEMS 

Since the onset of convective motions is closely related 
to the stability of the systems, it is worthwhile to consider 
the criterion on their stability. A general method for study- 
ing hydromagnetic stability problems has been established with 
respect to the non-rotating systems (Bernstein et al., 1958). 
The criterion on the onset of instability in such systems is 
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given by a following inequality: 


f ?* F(f) dv > 0 

In fact, any system which fulfilles this inequality is 

always unstable (Bernstein et al., 1958). As will be 

later discussed, § expresses a given displacement vector 

as applied to the system under consideration and is its 

— * — + 

complex conjugate. F(f) is given by equation (2-20) in the 

paper of Bernstein et al. (1958) and gives a generalized 

— > 

force induced by § (see (5) later). 


3. STABILITY IN THE ROTATING SYSTEMS 


When we investigate the stability of rotating hydrogmagnetic 
systems, we need consider the effect of the Coriolis force. 

Thus the discussion on stability of these systems is different 
from the forgoing one as made by Bernstein et al. (1958). 

For the rotating systems, the equation of motion is 
given as 


and 


P ^-| + 2p n x || = F (§) 
3 1 


(4) 


F(f)= grad [y P div § + (S’grad) P] + j x Q 


- B x curl Q + [div (p§)] grad 


( 5 ) 


where 5, t, f6 and y are the displacement vector due to 
perturbation, the time, the gravitational potential and 
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( 6 ) 


the ratio of two specific heats. Here Q(?) is given by 

Q(?) = curl (I x B) 

In the above equation, gradjrf consists of two terms 

as the gravitational and the centrifugal forces, but, in 

general, the contribution of centrifugal force is negligibly 

small (e.g., Chandrasekhar, 1961). In case which the system 

is rotating so fast like pulsars, the force F(§) necessarily 

2 

becomes a function of Q as suggested by Steinitz (1965). 

In this paper, we, however, will not consider such cases. 

— ► 

Here we assume that § is expressed as 

1 iwt „ v 

§ = § o (r) e , (7) 

’ ^ * 4 

where (r) is a complex positional vector. By substituting 

(7) into (4), we obtain 

- u) 2 p§ + 2iu) p (Q x I) = F(l) (8) 

By multiplying scalarly (a complex conjugate of §) with 

(8) , we futher obtain 

- 'JJ 2 pi*i + 2iuup i* (Q X i) = i*-F(|) (9) 

It is known that F(|) is self-adjoint and therefore the 
right side of (10) is always real (e.g., Bernstein et al., 
1958) . This equation is here integrated over the volume V, 
in which the system under consideration is involved, with the 
boundary condition as 

n-1 - 0, 
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where n is the unit vector perpendicular to the surface 
which encloses the above volume. Thus we obtain 

- uu^y*p §* § dv + 2oo y*p §* i (Q x §) dv 

=/ 1*-F(f) dv (10) 

It is clear that the first term on the left side of the 
above equation is real and positive. 

Since the vectorial operator i[ft x • •••] is Hermitian 
(see Appendix 1), the second term on the left side of (10) 
is also real. Thus, we can solve the above equation algeb- 
raically with respect to tu . 

The solution of (10) for id is given by 


p §*§ dv (*' 


i p §* (ft x |) dv 


"7 

+ ([ /i p S* (ft x §) dv] 2 - [/p in dv] [ F (§) dv] 


( 11 ) 


In order that id is real, it is necessary that 



F (I) dv ^ 0 


( 12 ) 


because 

/p dv > 0 and [ fi p §* (Q x §) dv] 2 > 0 

Inequality (12) gives a condition that the system under con- 
sideration is stable. In so far as this inequality is fulfilled, 
the system is always stable. 
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Even if inequality (12) is not fulfilled, the system 
is also stable when 

[ J' 1 p 1*(H x I) dv] 2 - [/p 1*? dv] [ / S*-F(I) dv] > 0. 

(13) 

In consequence, the system is unstable when the above two 
inequalities, (12) and (13), are not fulfilled: when in- 

equaltiy 

[ f±p l*(fi X I) dv yifp 1*1 dv] < dv (14) 

is fulfilled, the system is always unstable. The left side 

of the above inequality is always positive unless Q = 0. 

Although the form of § is not specified in the above 

inequality, the result of (14) gives a general criterion on 

instability of the system. Only when this inequality is 

fulfilled, the system is unstable. It should be remarked 

that the lowest value of the onset of instability increases 

2 

in proportion with Q as is deduced from the left side of 
(14). This means that the effect of the Coriolis force is to 
inhibit the onset of convection. 

4 . DISCUSSION ON THE ONSET OF CONVECTIVE MOTION 

Since the left side of (14) is always positive, it is 
clear that the effect of the Coriolis force is to stabilize 
the system by inhibiting the onset of convective motion. 
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In other words, the rotating system would, in general, 
be stabilized due to the action of the Coriolis force, 
because the lowest value for the onset of instability is 
raised by the amount of l J" ipS* (Ox §) dv]^/j- J* p§*"| dv ] 

(> 0) in comparison with that for the non-rotating systems. 

The stabilizing effect by the Coriolis force has also been 
found in the case for the rotating fluid heated below 
(Chandrasekhar, 1953). 

Once inequality (14) is fulfilled, the system moves to 
an over-stable state as is deduced from (11) . Then con- 
vective motion is set up in the system, but is oscillatory 
in nature. 

In these above discussions, we had better remarked 
some restriction on the applicability of the formulae con- 
sidered in the last section, because the present discussion 
cannot be applied to the case, in which the effect of the 
Coriolis force is so strong that the same perturbation § 
cannot be used for both the rotating and non-rotating systems. 
Namely, our present result is only useful for the rotating 
systems where the term of the Coriolis force can be treated 
as a small perturbation with respect to the non-rotating 
systems . 

As shown above, the onset of over-stability is usually 
inhibited by the effect of the Coriolis force. However, 
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once the state of over-stability is reached, convective 
motion is generated and, later on, seems to become turbulent 
due to the further growth of over-stability in the system. 
Since the build-up of convective motion is highly inhibited 
in comparison with the case for the non-rotating systems, 
it can be said that the rotating stars, in which there 
exists convective motion, are relatively rare. 

At present, we have known that, although stars are 
more or less rotating, all of them do not necessarily have 
their own magnetic fields. In the very early history of 
their evolution, the most stars might have captured the 
magnetic fields in the interstellar space and therefore 
they must have had their own magnetic fields during their 
young ages. As discussed earlier in this paper, these fields, 
however, might have decayed throughout their evolution by 
the Joule dissipation unless there existed convective motion 
in their interior. Since the effect of the Coriolis force 
is to inhibit the onset of convective motion, the rotation 
of stars must have worked as an obstracle on the maintenance 
of their own magnetic fields while the systems are stable. 

This fact seems to be one of the reasons why, at present, 
very limited number of stars and planets only have their own 
magnetic fields. 
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For the first time, the effect of the Coriolis force 
is to inhibit the onset of convective motion. Once this 
motion- is set up, it seems, however, that the role of this 
force enforces to build up turbulent motion as a result of 
over-stability. Important is that turbulent motion in the 
rotating systems is non-re flectionsymmetric on account of 
the action of this force (e.g., Braginskii, 1965; Steenbeck 
and Krause, 1969). Thus the onset of turbulent motion in the 
rotating systems seems to be very important from the view 
point on the dynamo action for the maintenance of the magnetic 
fields of stellar objects (e.g., Parker, 1969, 1970). 

5. CONC LUDING REMARKS 

In this paper, we have shown that the effect of the 

Coriolis force is to inhibit the onset of convection in the 

rotating hydromagnet ic systems. Since, as is evident from 

(14), the lowest state on the onset of convective motion is 

2 

proportional to Q , the effect of the Coriolis force becomes 
larger as the magnitude of Q increases. In this case, the 
state which is not stable is always over-stable. Turbulent 
motion which seems to develop from the oscillatory growth 
of convective motion would be effective on the dynamo action 
on magnetic fields in the systems. 
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APPENDIX 1 


The second term on the left-hand side of (10) is re- 


written as 


fl* i 


i(Q x §) dv 




§* i(e x §) dv, 


(A-l) 


where e is the unit vector in the direction of Cl. 


Let us examine whether the above term fulflles the follow 


ing equality 


^ cp A i/ dv =y A ijr cp dv, (A-2) 

where cp and \|r are functions of complex variables and A is 
the operator (see Schiff, 1968). In our case, 


cp = I*, i i = i 


and 


A = i [e x* • • 


By using these definitions, we calculate the right-hand 
side of (A-2): 


- 0/i(e x I) * I dv 

~ Cl f 1 I (ex §) dv = Cl /i(e x §) § dv 


f 1 


‘f 


Ci I i (ex §) dv 


(A-3) 


This calculation thus shows that expression (A-l) is 
Hermitian: i.e., the operator A is Hermitian. This means 

f - Z - . 

that Cl I (e x §) I dv is purely imaginary. 
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